Harmonic lattice behavior of two-dimensional colloidal crystals 
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Using positional data from video-microscopy and applying the equipartition theorem for harmonic 
Hamiltonians, we determine the wave-vector-dependent normal mode spring constants of a two- 
dimensional colloidal model crystal and compare the measured band-structure to predictions of the 
harmonic lattice theory. We find good agreement for both the transversal and the longitudinal mode. 
For q — > 0, the measured spring constants are consistent with the elastic moduli of the crystal. 

PACS numbers: 61.66.-f,63.20.Dj,82.70.Dd 



Colloidal crystals have a long tradition as condensed 
matter analogs of ordinary solids. They are studied, 
for example, to understand phenomena such as freezing 
and melting [lj. Unlike in ordinary solids having prop- 
erties that are often difficult to connect to the under- 
lying atomic interactions, the interparticle potentials in 
such colloidal crystals are in most cases precisely known 
and, more importantly, externally controllable. More- 
over, the relevant time- and length scales in colloidal 
systems are comparatively easy to access experimentally. 
Both aspects suggest studies directly probing the con- 
nection between microscopic interaction-potentials and 
macroscopic crystal properties. 

The property we here consider is the crystal's elastic 
response to thermal excitations, specifically, the phonon 
dispersion relations. In this regard, colloidal crystals are 
rather special in that their phonons are almost always 
overdamped: the ratio between the wave-vector depen- 
dent frequency ui(q) = y/X{q)/m, characteristic of the 
harmonic forces with spring constants X(q), and the fric- 
tion factors A(q) (also g-dependent 0) for the modes 
of lattice motion through the host liquid, is typically 
of the order 1CU 3 to 1CU 4 in colloidal systems. There- 
fore, the time autocorrelation function of a phonon nor- 
mal mode coordinate decays exponentially with a rate 
given by X(q)/A(q) 0,13- This decay rate can be, and 
has repeatedly been, measured by means of dynamical 
light scattering 0, 0, or inelastic light scattering 
[8j . Phonon dispersion relations have been determined in 
charge-stabilized 0- 0- E3| and purely hard sphere colloidal 
crystals 0, , in the context of dusty plasma physics , 
but also in more exotic systems such as crystals made of 
mm steel spheres or optically anisotropic spheres . 

Microscopic information about the spring constants 
and thus the particle interaction potentials can only be 
derived from these decay rates, i.e. from X(q)/A(q), if 
one resorts to a model describing the complicated fric- 
tional forces, especially those of hydrodynamic origin. A 
direct access to X(q), i.e. one free from any assumptions 
of a model, is not possible in this approach. 

In this Letter we report on a video- microscopy study of 
two-dimensional (2D) colloidal crystals and show how to 
obtain direct access to the normal mode band-structure 



FIG. 1: Micrograph (420 x 310 /im) of a typical colloidal 
crystal investigated in the present work; the two-dimensional 
system consists of paramagnetic colloids confined at the air- 
water interface of a hanging water drop. 



X(q) of the crystal, circumventing, in particular, the diffi- 
culties arising from the hydrodynamic interactions. The 
central idea is to avoid a dynamical measurement and 
to analyze instead spatial correlations between the par- 
ticles which are then related to the X(q), the eigenvalues 
of the dynamical matrix characterizing the elastic prop- 
erties of the harmonic crystal. This becomes possible 
through use of digital video-microscopy [TH providing us 
with the trajectories of all particles. The colloidal system 
we examine is well-studied and the interparticle potential 
precisely known [l2l Hij , . This will be of advantage 
when establishing a quantitative link between the mea- 
sured X(q) and the theoretical band-structure based on 
the pair-potential. 

The experimental setup is an improved version of 
the one described in [l^: Spherical colloids (diameter 
d = 4.5 /im) are confined by gravity to a water/air inter- 
face formed by a water drop suspended by surface tension 
in a top sealed cylindrical hole of a glass plate. The flat- 
ness of the interface can be controlled within +/- half 
a micron. The field of view has a size of 835 x 620 /im 
containing typically up to 3 • 10 3 particles, whereas the 
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whole sample contains about 10 5 particles. The num- 
ber of particles in the field of view is controlled by the 
curvature of the droplet via an active regulation with an 
accuracy of 1% and the biggest observed particle-density 
gradient is 0.7%. The variation of the inclination of the 
sample is in the range of a « 1 /iRad so that the col- 
lective motion of the particles relative to the field of 
view is below 2 pm/h providing best equilibrium con- 
ditions for long time stability. The particles are super- 
paramagnetic due to Fe2C>3 doping. A magnetic field 
B applied perpendicular to the air/water interface in- 
duces in each particle a magnetic moment M = X-B 
which leads to a repulsive dipole-dipole pair-interaction 
energy of (3v(r) = T/( A /irpr) 3 with the dimensionless in- 
teraction strength given by T = /3(/io / '4tt)(x-B) 2 p) 3 ^ 2 
([3 = 1/kT the inverse temperature, \ the susceptibil- 
ity). The interaction can be externally controlled by 
means of the magnetic field B. T was determined as 
in Ref. (0|) an d is the only parameter controlling the 
phase-behavior of the s yst em. For V > 60 the sample 
is a hexagonal crystal [13j (see Fig. 0J). The sample 
was tempered at high interaction strength up to T = 250 
deep in the crystalline phase until a 2D-mono-crystal was 
observed. We here analyze three different crystals, from 
hard to soft (r = 250, 175, 75), and use for each system 
about 2000 statistically independent configurations with 
approximately 1300 particles, recorded at equal time in- 
tervals (At = 2s) in a 440 x 440 fim frame using digital 
video-microscopy with subsequent image-processing on 
the computer. For each of all N particles in a given con- 
figuration, we determine the displacement u(R) of the 
particle from its equilibrium position R. 

Using the theory of harmonic crystals ^5|, we now 
derive an equation guiding us from the measured dis- 
placement vectors u(R) to the eigenvalues of the dynam- 
ical matrix. Let D^^q) (/i, v £ {x,y}) be the dynam- 
ical matrix |l5| . connected through a Fourier transfor- 
mation to the matrix D^^(R,R') which is essentially 
the matrix of the second derivatives of the pair-potential 
v(r) ~ r/r 3 . It is obvious that D^^(q) depends linearly 
on the interaction strength parameter T; therefore we 
write D^^ v (q) = (kTT /a 2 )D fliV (q) and obtain the dimen- 
sionless dynamical matrix D^^q) which is independent 
of r (a is the lattice constant of a hexagonal lattice). 
From the measured u(R) we seek to determine the eigen- 
values of D^ v {q) which we denote by X s (q)a 2 /(kT T). 
Here the polarization subscript s stands for the longitu- 
dinal (s = I) and transversal (s = t) mode. The harmonic 
potential energy of the crystal can be written in the fol- 
lowing form 
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placement vectors u(R). The equipartition theorem 
for classical harmonic Hamiltonians states that on av- 
erage every mode has an energy of kT/2. Thus 
{\/2V)(ul{q)D il , u (q)u u (q)) = kT/2 and this leads us to 



±(u;(q)u„(q))=kTD-i(q) 



(2) 



where in our case the average has to be taken over all 
measured configurations. Introducing with p s (q) (s — 
i, I) an abbreviation for the eigenvalues of the matrix 
T(u* (q)u v (q))/ (Va 2 ), one arrives at 
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with V = Nv = Ny/3a 2 /2 and u v (q) being the 
f'th component of the Fourier transform of the dis- 



(s = t,l). 

Static and slowly moving distortions of the lattice are 
the main source of error in our experiment. We some- 
times observe long-range bending of lattice lines, a finite- 
size problem which in soft crystals can be partly overcome 
in giving the crystal enough time to equilibrate. This 
takes the more time the harder the crystal is. For our 
hardest crystal we have not managed to avoid a small, 
but clearly visible bending of lines. A second problem 
is related to the determination of each particle's equilib- 
rium position R, without which the displacement vectors 
u(R) cannot be determined. A cooperative drift of all 
particles can be observed, a behavior typical of 2D crys- 
tals. An illustrative example pictures can be found in |l7| 
and also in 0|: particles can depart significantly from 
their lattice sites, but keep a nearly constant distance 
from each other so that nearby trajectories are similar. 
However, when resorting to a local coordinate system 01 
the overall symmetry remains crystal-like, and the crys- 
tal melts in accordance with the Lindeman criterion in 
spite of the collective drift (measured in the global coor- 
dinate system, the root mean square displacement would 
diverge at long times, indicating the instability of the 2D 
crystal \La[). These observations have been confirmed ex- 
perimentally [l3| . To correct our data for drift, we first 
calculated coarse-grained trajectories by averaging over 
a sliding time window AT having a width of 25, 40 and 
60 At for T = 250, 175 and 75, respectively. We then an- 
alyzed the short-time displacement of the particles with 
respect to these coarse-grained trajectories to obtain the 
true fluctuations of the underlying crystal. Our whole 
data evaluation procedure has been successfully tested 
by processing data obtained from Monte Carlo simula- 
tions, using the pair-potential and parameters of our ex- 
periment. 

Fig. (0) shows l/p s (q) from eq. 10 as obtained from 
the measured set of displacement vectors for V — 250, 
175 and 75, and compares it to the theoretical band- 
structure (solid lines) of a harmonic crystal having a 
two-dimensional hexagonal lattice (a = 12.98 fim). The 



FIG. 2: Band-structure of harmonic lattice spring constants 
of a 2D colloidal crystal. Symbols for constants experimen- 
tally determined with eq. and from the relative dis- 
placement of the particles from their equilibrium position 
for a soft (F = 75, empty circles), hard (F = 175, empty 
squares) and very hard (F = 250, filled triangles) crystal; 
solid lines for the theoretical band-structure calculated from 
standard harmonic crystal theory using the pair-potential 
f3v(r) — Y /{y/Wpr) 3 . Fhe inset shows the first Brillouin zone 
of the hexagonal lattice and labels for high-symmetry points, 
defining the lines in the interior and on the surface of the first 
Brillouin zone along which the band-structure is plotted. The 
upper curve corresponds to the longitudinal, the lower one to 
the transversal mode. 



latter is based on the second derivatives of the known 
pair-potential and results from diagonalizing D^ u (q) 
[15|. 17 neighbor shells have been taken into account 
in D fJ _, u (R, R'), the difference to the results for only 3 
shells is already tiny. We find good agreement for both 
the transversal and longitudinal mode. No fit param- 
eter has been used. Pre-averaging with a finite time- 
window improves the agreement. Without it (i.e. if tak- 
ing the average over the whole trajectory to define R), 
the peak at the M-point in the band-structure is about 
10 % smaller for each T than it is in Fig. (0. The data 
are particularly sensitive to the quality of the crystal near 
the edges of the first Brillouin zone, especially near the 
Af-point. The uncertainty in determining the direction 
of the lattice lines plus the bending of these lines explain 
the remaining differences between the theoretical and the 
experimental band-structure, but also the differences be- 
tween the three different crystals. We also checked for 
the occurrence of dislocations in all our samples. Only 
thermally activated dislocation pairs have been observed, 
but no static, isolated dislocation destroying the crystal 
symmetry. 

It is instructive to study the q — ► limit. We define 
q\\ = q x e x + q y e y and q± = -q y e x + q x e y where q = 
Qx&X "T qy e y , and multiply on either side of eq. @ q\\ and 



FIG. 3: d||(g)r and d±(q)T from eq. Q, evaluated with 
the same data as taken in Fig. @. For q — > 0, d\\(q) — > 
kT/((K + n)a 2 ) and d±(q) -> kT/(fia 2 ) where K is the bulk 
modulus and n the shear modulus. Curves (1) and (2) are 
d±(q)T in the r -> M (curve (1)) and the F -> K (curve (2)) 
direction, (3) and (4) are d\\(q)/T in the F — ► K (curve (3)) 
and the F — > M (curve (4)) direction. Symbols, lines and the 
labels F, K, M for symmetry-points are all defined in Fig. (|5J. 
Arrows indicate the T = prediction of the elastic moduli, 
see text. 



gj_ from the left and the right, to find 

di(q) = y(\qiu(q)\ 2 ) = kT^q^D^^q^u, i = \\,±. 

(4) 

It can be shown that lim g ^ | = kT/((K + /u)a 2 ) 
and \im q _> d±(q) = fcT/(/ia 2 ) where K and /i are 
the bulk and shear elastic moduli of continuum the- 
ory. Since D^ M (q) depends linearly on T, the quantity 
Tdi(q) (i = ||, _L) is independent of T. Fig. @ shows 
<i|| (q)T and d±(q)T, evaluated with the experimental data 
for all three crystals (T — 75, 175, 250, symbols) and 
theoretically using the dynamical matrix D^^(q) and 
the pair-potential (solid- lines). For the pair-potential 
~ T/r 3 the elastic constants can be calculated to be 
Ka 2 /kT = 3.46ir and fj, = K/10 in the limit r — > oo 
(T = 0) [3 G3. Arrows in Fig. © indicate the pre- 
diction of the T = calculation. While the data on 
the longitudinal branch (curve (3) and (4)) show excel- 
lent agreement and correctly approach the T = inverse 
bulk modulus, theoretical and experimental data on the 
transversal branch disagree at low q. This is clearly a 
finite size effect. Taking data from this very experiment, 
it has already been shown 14] that appropriate finite 
size scaling leads to an almost perfect agreement with 
the T = 0-prediction of the elastic moduli. Wc should 
also remark that the location of the branching points of 
the r — > K and r — > M curves in Fig. ® reveal that 
the assumption of isotropy is justified only if qa > 1 for 
the transversal and qa > 3 for the longitudinal modes. 
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This defines the limits for the continuum approach often 
chosen to describe this system. 

In summary, we used video-microcopy data to deter- 
mine the wave-vector dependent normal mode spring 
constants of a two-dimensional colloidal model crystal. 
We checked the continuum limit and compared the 
experimental data to the predictions of the classical 
theory of a harmonic crystal. Our data evaluation 
procedure can be seen as an illustration of the validity 
of the equipartition theorem which we used to derive 
the spring constants from the particles' displacement 
vectors. The success of our undertaking was not clear 
from the beginning; if a crystal in 2D is not stable, how 
can one measure the normal mode spring constants? 
Here the ideas put forward in 0, 0] proved help- 
ful, specifically the introduction of a local coordinate 
system. Analyzing particle distributions at equilib- 
rium, we were allowed to completely ignore the lattice 
dynamics. In this context, it is worth remembering 
that Fig. (J2J) is not a phonon-dispersion relation in the 
classical sense as there are no phonons propagating 
with u> s (q) — \f \ s {q)/m. Our results suggest that it 
is more appropriate to think of a colloidal crystal as a 
bead-spring lattice immersed in a viscous fluid ji^. A 
normal vibration mode then transforms into a 'normal 
relaxation mode' 0, Ei3- an d the motion of a particle 
is to be understood as superposition of these 'normal 
relaxation modes'. A time-dependent analysis of our 
data would allow to study the relaxation process of these 
normal modes. From the relaxation times and armed 
with the spring constants measured here, one could then 
proceed to evaluate the (f-dependent friction factors to 
study, for example, hydrodynamic forces. These avenues 
await further investigations. We finally remark that 
the statics and dynamics of overdamped phonons in 
two-dimensional colloidal crystals may also be seen as 
a contribution to our understanding of surface phonons 

m- 
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